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ABSTRACT 

Galaxy merger simulations have explored the behaviour of gas within the galactic 
disk, yet the dynamics of hot gas within the galaxy halo has been neglected. We report 
on the results of high-resolution hydrodynamic simulations of colliding galaxies with 
hot halo gas. We explore a range of mass ratios, gas fractions and orbital configurations 
to constrain the shocks and gas dynamics within the progenitor haloes. We find that : 
(a) A strong shock is produced in the galaxy haloes before the first passage, increasing 
the temperature of the gas by almost an order of magnitude to T ~ 10^-^ K. (b) The 
X-ray luminosity of the shock is strongly dependent on the gas fraction; it is > 10^^ 
erg/s for halo gas fractions larger than 10%. (c) The hot diffuse gas in the simulation 
produces X-ray luminosities as large as 10^^ erg/s. This contributes to the total X- 
ray background in the Universe, (d) We find an analytic fit to the maximum X-ray 
luminosity of the shock as a function of merger parameters. This fit can be used in 
semi-analytic recipes of galaxy formation to estimate the total X-ray emission from 
shocks in merging galaxies, (e) ^ 10-20% of the initial gas mass is unbound from the 
galaxies for equal-mass mergers, while 3 — 5% of the gas mass is released for the 3:1 
and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible 
mechanism to enrich the IGM with metals. 
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1 INTRODUCTION 

Galaxy mergers drive galaxy evolution; they are one mechanism by which galaxies ^row. The merger process is also accompa- 
nied by a strong gas inflow towards the central re^o n ( Barnes &: Hernquistll 1991 : Mihos Hernquistll 19941 : Combes et al.|[l994 : 
Hernquist k Mihosl[l995l : ISanders k Mirabell[l996l ) triggering a central starburst and perhaps an AGN. Consequently, galaxy 
mergers have been studied quite extensively both numerical l y (lHolmberdll 941: Toomre & T oomrelll972l : iBarnes k Hernquist! 
1992. ll996l:lMihos Hernouistll996l:ISDringel Whitelll999l:lDubinski et ah 1999: Barnes 20oi l and observationallv (|Schweizer 
1986 : lRoberts et all 20021 : Laine et all 20031 : Casasola et al. 20041 ) . However, none of thes e galaxy merger simulations included 
the h ot gas in the galactic halo that is expected from standard galaxy formation theory ( White ReesI 19781 : White Frenk 
199ll ). In this paper, we explore the effect of hot halo gas during galaxy mergers with high-resolution numerical simulations. 

Semi-analy tic galaxy forraation models have long pre dicted the existence of an extended reservoir of hot halo gas around 
spiral galaxies ( White k ReesI 1978 : White Frenk|[l991 ). In this scenario, gas falling into the dark matter potential shock- 
heats to the virial temperature, Tvir = lO^K (vcirc/167km/s)^ and subsequently cools over a characteristic time-scale, tc — 
kT/niA(T), that depends on the temperature T, the number density rii and the cooling function A(T). This gas cools via 
thermal bremsstrahlung and atomic line emission processes, with the majority of the radiation occurring in the soft X-ray 
band. The expected X-ray luminosity, Lx , is a steep function of the circular velocity, Vdvc - Lx oc v^irc- Since the effect of 
mass (circular velocity) on this halo X-ray luminosity is so non-linear, any attempts to observe the X-ray signature of this gas 
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would naturally be biased towards very massive haloes. Consequently, hot ^as in massive elliptical gal axies has been shown by 
X-ray observations of cluster galaxies as well as in isolated ellipticals (see Mathews Brighentil 20031. and references therein). 
Observations of massive ellipticals show that the X-ray luminosity is non-stellar in nature (|0 'Sullivan et al" 200 ih leading 
credence to its hot gaseous or iffln . An extended gaseous halo in a quiescent spiral has been observed only for the massive spiral 
NGC 5746 JPedersen et al.ll2006l ll% owever, the existence of a hot, rarefied an d extended galactic h alo has been inferred from 
OVI absorption lines in high veloc ity clouds in the vicinity of the Milky Way (|Sembach et al.ll2003l ). The density of this halo 
gas (riH ~ 8 X 10~ cm" (M urahll200d : ISnowden et al.lll997l ) is consistent with the predicted density for a Milky Way type 
galaxy. Given the existence of hot halo gas, it is important to include it in galaxy merger simulations. 

During a merger, the bulk motion of the galaxy gas is converted into thermal energy via shocks. The orbital velocity of 
the gas during a merger is ^ 100 km/s whereas the sound speed is ~ 10 km/s. This can produce strong shocks with Mach 
numbers > 10. Assuming an adiabatic index of 5/3, the shocks may heat the gas to: 
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where Tfinai is the gas temperature after the shock, /j is the mean molecular weight, mn is the proton mass, k is the Boltzmann 
constant, and Vinitiai is approximately equal to the radial velocity prior to the collision. Nearly 50% of the kinetic energy can 
be converted into thermal energy of the gas particles. Thus, we expect that the hot gas will shock-heat and radiate in soft 
X-rays. The shocks may occur as soon as the gaseous haloes collide - well before the optical extents of the galaxies begin to 
show signs of morphological disturbance. 

This hot halo gas may play an important role in enriching the intergalactic medium (IGM). Galaxy formation theory 
predicts that baryons fall into the dark m atter potential and shock-heat. Absorption l i nes reveal th e presence of diffuse 
gas with column density A^h ^ 10^^ /cm^ 



ransson EDsteinlll982l: ICowie et al.lll995l : IWand Il995h . At z ^ 3, most of 



the baryonic mass is in the Lyman-a forest ( Fukugita. Hogan. PeeblesI 1998 ). The baryon fraction in the Lyman-a forest 
systematically decreases with z and is almost zero at present. It would be reasonable to expect that all the baryons are 
now part of galaxies in the form of stars and the interstellar medium. However, recent estimates of the massfraction of 
baryons in different comp onents still fail to d etect ^60% of the total baryon content at low with 6% i n stars in ga laxies 
( Fukugita Peeblesll2004l ) . 2% in neutral gas ( Zwaan et al.l 2003), 4% in galaxy clusters and galaxy groups ( Allen et al.| [2002) 
and 29 ± 4% in \ow-z Lyman-a absorbers (|Penton et al.ll2004l : IPanforth Shulll l I2OO8I ) . The consensus is that the majority 



of the 'missing' gas is located in regions of low over-density, Sp/p ^ 10° — 10^ with temperatures in the range 10^ - 10^ K. 
This is the Warr n-Hot Intergala ctic Medium (WHIM) and is currently one of the most elusive astrophysical objects (for a 
recent review see iBregmanI 120071). Numerical simulations pre d ict that the WHIM must reside in filaments and r elatively low 



over- density regions (ICen Ostrikei]ll999l : IPave et al.l[l999l : ICen Ostri'tol I2OO6I : iDave Oppenheimerl I2OO7I ) . Numerous 



efforts have been taken to detect and characterise the WHIM; however, most of the observations to date are not iron-clad and 



have required some fine-t uning in the fitting procedure for the data to come up with a significant detection (|Nicastro et al. 



2OO5I : Mathur et al. 20031 ) . The current instrumentation limits on the Chandra and XMM-Newton are a major hold-up in the 
discovery of the WHIM; next generation X-ray telescopes like Constellation-X should be capable of detecting the absorption 
lines from the WHIM along the line of sight to bright background sources. 

The immediate question is: How does gas get into the WHIM state ? Some form of mass and energy injection is essential 
to create this hot reservoir of gas, and this form of feedback must also account for the metals in the IG M. Numerical simu- 
lations of structure formation must incorporate various baryonic processes such as supernova feedba ck (|Mac Low Ferrara 
I999I : iBorgani et alll2002l: Irha cker Couchmanl'200d: 'Scannauie co et al.ll2006l : iGovernato et al.ll2007l) , starburst-driven galac- 



tic winds (ISpringel Hernau ist 2003; Oppenheimer & Dave 200 8|) and AGN feedback (ISpringel et al, 



2OO5I : IScannapieco et allboosl : iThacker et al.ll2006l : iHopkins et al.ll2006l : ICroton et al.ll2006l : ISiiacki et al 



I2OO5I : iDi Matteo et al 



20071: iBhattacharva et al 



2008 : 151 Matteo et ahlbo osl) . All these energy delivery mechanis ms are more effective in low mass galaxies ( IScanna pieco et al 
2008), e.g., these prescriptions blow-out the gas in dwarf galaxies ( Mac Low &: Ferrarall 19991 : Kobavashi et a l. 2007) . Howe ver, 
this is in direct disagreement with observational results that show dwarf galaxies to be more gas- rich (|Kannappanll2004l ). In 
this paper, we will explore if galaxy mergers, which are ubiquitous, can be used to fuel the IGM. 

Observation s of major galaxy mer gers show that the merger rate evolves rapidly, increasing from ~ 5-10% at ^ 1 to 
~ 50% ai z ^ 3 (|Conselice et al.ll2003l ). Integrating this merger rate suggests that a typical massive galaxy will undergo 4-5 
major mergers between ^ ~ 3 and z ^ 0, with most of these mergers occurring at z > 1.5. Nuni erical simulations can be used 
to explore the evolution of the merger rate as well as the merger ratio with z. IPakhouri Mai (2008) analyse the Millenium 
Simulation to constrain this evolution and find that a major merger (merger ratio ^1/3) occurs at the rate of 0.5 per halo 
per unit redshift. Minor mergers with ^ ^ 1/10 are far more common at all redshifts. 

In this paper, we explore typical mergers with a Milky Way type primary galaxy. We sample the major mergers by 



^ It has since been shown that the detection of hot gas in NGC 5746 may be caused by contamination from unresolved point sources (J. 
Rasmussen, private communication). 
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simulating an equal-mass merger between two such galaxies. We also attempt to get the boundary of the major mergers for 
a merger ratio, ^ 1/3 for our second set of merger ratios. Minor mergers are the most frequently occurring mergers in the 
Universe; we model this by simulating a merger with ^ ^ 1/10. This describes the three independent sets of simulations that 
we will perform to constrain the effect of the merger ratio on the behaviour of the hot gas during galaxy mergers. Mergers 
need an additional specification - the merger orbit. We will describe the parameters for specifying the content of the galaxies 
(see Sections 12.11 and |2.2|) . and construction of the merger orbits in Section \2A\ In Section |4] we will present the results of 
our numerical simulations of galaxy mergers. We will also derive an analytic fit for the peak shock X-ray luminosity and 
the unbound gas mass as a function of the progenitor galaxy masses and impact parameters. In Section [5] we discuss the 
implications of our results and future work. 



2 METHOD 



We create an equilibrium model for an isolated galaxy by numerically sampling for the joint distribution function (hereafter, 
DF) of the dark matter and gas. The DF for a spherically symmetric potential-density pair can be written as: 



/(f) = 



V87r2 



/' 

^0 



(2) 



d^2 y^TT^ 

where, 8 and ^ are the absolute value of the total energy and the gravitational potential respectively and p is the density. 
We will describe the initialisation of the dark matter and the gas component in the following subsections. 



2.1 Constructing a dark matter halo 



Large scale cosmological simulations lead to a two-parameter universal dark matter halo density profile as shown by Navarro, F renk. &: White 
(fl997) (hereafter, NFW). The key characteristics are the virial radius, Rvir such that the average density of the halo within 
that radius is a certain multiple (A200 = 200) of the critical density of the universe, pcrit = SHq/SttG, and the concentration 
parameter, c = Rvir/rs, where rs is the scale radius of the NFW profile. Since the mass of the NFW halo is infinite, we 
implement an ad- hoc truncation radius for the dark matter (and the gas) at i?haio = 1.2 x Rvir- 

To compute the DF, we first create a logarithmically spaced array in r with 10^ bins. The minima and maxima of the 
bins are fixed at 10~^ and 100 times the virial radius. Thus, we capture the inner cusp of the NFW halo profiles as well as the 
outlying particles with highly eccentric orbits; this can be seen from the upturn in the DF values for small E in Fig. [T] Note 
that the inner limit is much smaller than our softening length while the outer limit is much larger than the truncation radius 
of 1.2. This was done to ensure a higher accuracy of the DF itself, even though at present we can not reach such numerical 
resolution. On this finely-spaced grid, we define the values for p(r), M(r) and $(r) and obtain the derivatives by finite 
differenci ng. The DF is then obta ined by numerically integrating Eqn.[2]with a using a fourth-order Runge-Kutta integrator. 
Following Lokas & MamonI (|200ll ). we chose our system of units as C = 1.0, Mvir = 1.0, Rvir = 2.0 with the resulting f{8) in 

units of VSMvir/ (Rvir Vvir)^. 

Our fiducial galaxy is a Milky- Way type galaxy with a circular velocity of 160 km/s and a concentration parameter c = 10. 
Mvir for this primary galaxy is 1.36 x 10^^ M© with Rvir = 228.6 kpc. The secondary galaxies for our merger simulations are 
chosen to be one-third and one-tenth of the mass of the primary; this translates approximately to c ircular velocit i es of 1 10 and 
74 km/s. To determine the concentrations of such haloes, we use the analytic fitting formula of iBullock et al.l (|200ll ). With 
this prescription, our fiducial galaxy has a concentration of c = 10, and the concentrations for the two secondary galaxies are 
16 and 25 respectively. 

Fig. [1] shows the comparison between the DF from ILokas MamonI ( 200lh for an NFW halo with c = 10 and the one 
obtained with our code for the fiducial galaxy with a 10% gas fraction. 



2.2 Constructing a gas-rich galaxy halo 



Now that the dark matter has been initialised, we need to initialise the gas. Earlier papers have used r^b/^^dmPdm as the 
ad-hoc density for the gas for cosmological simulations of structure formation. While that approximation is valid for initial 
conditions in a cosmological framework where the force field is completely determined by the dark matter (and gas merely 
falls into the dark matter potential), our simulations are in a fixed s pace-time and we needed an analy t ic physically motivated 
gas density pro file. We chose the observationally verified /3-profile ( Cavaliere Fusco-Femianol 19761 : Jones Forman 1984 : 



Eke et al 



19981 ) used to model gas in galaxy clusters, a more massive analog to our gaseous haloes. The gas is isothermal and 
in hydrostatic equilibrium. We use /3 = 2/3 and a core of rc = Ts/3, giving a density profile: 

Po 



Pgas 
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Figure 1. The c omp uted DF for our fiducial galaxy with c = 10 and 10% gas fraction. Overplotted in magenta is the DF from 
I L okas &: MamonI (|200 1) for a dark matter only NFW halo with Vvir = 160 km/sand c = 10. The presence of the gas makes very little 
difference in the overall DF. 



Since the density profile is spherically symmetric, we employ the Eddington inversion method to compute the DF for the joint 
dark matter and gas distribution. For the dark matter component, position- velocity pairs are drawn from this DF. For the 
gas, we chose the position from the gas density. However, the gas velocity depends on the temperature and the temperature 
itself can not be freely assigned since it is determined by the overall gas density-potential pair. We assume that the gas is 
polytropic and follows an equation of state, P oc p^^^ where P and p are the pressure and density of the gas. Finally, we 
assume hydrostatic equilibrium to compute the temperature of the gas self-consistently. The gas is assumed to be at virial 
temperature at the virial radius with a central density constrained by the total amount of gas in the halo. Once the temperature 
has been assigned to the gas particles, we calculate the velocity dispersion by assuming virialization. The final gas velocities 
are drawn from a Gaussian distribution with zero mean velocity and a velocity dispersion given by — 2kT/('y — l)pmH- 

The input parameters to our galaxy initialisation code are the circular velocity of the halo, the number of halo particles, 
the gas fraction and the number of gas particles. Our fiducial galaxy is similar to the Milky Way with Vvir = 160 km/s. 
To investigate the role of mass-ratio during mergers, we constructed galaxies with 1/3 and 1/10 of the mass of our fiducial 
galaxy. For a galaxy with one-third the mass, the virial velocity is 110 km/s and for one-tenth the mass, the virial velocity is 
74 km/s. Since we aim to simulate the mergers of two different mass galaxies, it was imperative to have the same mass per 
particle between each of the galaxies to prevent spurious numerical relaxation. Thus, the number of particles in the secondary 
galaxies are set by the corresponding number in the primary galaxy. We use a fixed number of 5 x 10^ particles for dark 
matter and gas in the primary galaxy for all the gas fractions. Since the total halo mass of the primary galaxy is fixed, our 
halo mass resolution is 3 x 10^ M©, independent of the gas fraction. The gas mass resolution for the 1%, 10% and 18% gas 
fractions in the primary galaxy is 1.5 x 10^, 1.5 x 10^ and 2.7 x 10^ M© respectively (see Table [T])0 

Since we are exploring how gas-rich galaxies behave during galaxy mergers, we did not want the input temperature 
profile, T(r), to vary and devised a multi-step process to reduce this temperature evolution. Initially, the gas temperature is 
determined from the hydrostatic equilibrium condition in the joint dark matter and gas potential. Then the gas velocity is 
assigned from the velocity dispersion based on virial equilibrium. The problem is that this assumption is not valid in the inner 
regions of the galaxy; even in isolation, the central gas particles shock-heat, increase their temperatures and decrease their 
velocity dispersion rapidly. A higher artificial viscosity parameter (AVP) aids this transformation and the core temperatures 
can increase by as much as a factor of 2 from the originally solved temperature profile. To alleviate this problem, after 
initialising the galaxy with the method described above, we instead evolve the galaxy model in isolation with AVP of 0.5 
to allow the gas particles to attain a new equilibrium temperature profile (see next section on the choice of AVP). This is 



^ Note that this gas mass resolution is the product total number of neighbours of a gas particle, 50, and the mass of each gas particle. 
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the first iteration on an equilibrium temperature profile for the halo gas. We find that the new temperature profile is indeed 
larger at the centre and turns over at the outer edge with an exponential drop-off near r ^ O.SRvir- We fit this temperature 
profile and apply it as a second iteration on the initial temperature profile for our model (see Fig. [2] for a comparison of 
the two iterations). This evolved temperature profile is nearly flat to ^ 0.5Rvir and is consistent with o bservations of galaxy 
clust e rs showing that the t emperature stays to within 40% of the central temperature out to 0.6 Rvir (|Leccardi &: Molendi 
20081 : IVikhlinin et al.1 liooil l. Ef inally, we evolve the galaxy with zero artificial viscosity over a period of ^ 1 Gyr to get 



rid of dynamical irregularities in our temperature profile fit (see Fig. [3]). Since there is absolutely no entropy increase with 
zero artificial viscosity, our temperature profile essentially remains constant in this step. We re-calculate Mvir, Rvir and Lx 
(see Table [2]) after iterating the initial conditions. For the primary galaxy, the maximum change in Rvir is ^ 2.5%; for the 
one-third galaxy and the one-tenth galaxy, the corresponding changes are ~ 3% and ^ 2.6% respectively. Once the galaxies 
are initialised, we check for stability by evolving them in isolation. One of the key ways of ascertaining the equilibrium of 
mult i- component systems is by computing the Virial of Clausius, VC. An equilibrium model for a galaxy should have VC 
fairly close to 1. The numerical models for our galaxies in the first iteration have a VC ~ 1.2, irrespective of the gas fraction 
or the halo size. Thus, while being close to equilibrium, the galaxies are not exactly in steady state initially. The value of VC 
for the isolated galaxies with the iteratively-obtained temperature profile, ranges from 1.12 for the primary to 1.07 for both 
the secondary galaxies. 

To test the stability of this new galaxy model after iterating the temperature profile, we simulate this evolved galaxy 
with AVP = 0.5 for a Hubble time. We find that the change in the overall temperature profile drops t o ^ 10 — 20% over 



this entire timespan. Our choice of softening length was guided by the suggested optimal softening from iPower et al. 
with Copt ^ 4 Rvir/VA^vir whcrc A/'vir is the number of particles within the virial radius. Our e = 1.1 kpc is comparable to the 
optimal value of Copt ~ 0.9 kpc. 

We use the parallel, hydrodynamic code GADGET-2 ( Springel et al. 2001; Springell 2005h for all the numerical sim- 
ulatio ns in this paper. GADGET-2 computes the gravitational forces with a hierarchical tree algorithm (! Barnes Hut 
1986h_ while ga s particles receiv e additional hydrodynamic acceleration as calculated using Smooth Particle Hydrodynam- 



ics (jOingold Monaghanl 1977 , hereafter SPH). The code explicitly conserves energy and entropy for the SPH particles and 
uses adaptive time-steps for the time-evolution. Since this is a first attempt to model the gas behaviour during the merger, 
cooling, star formation and radiative feedback have all been neglected. 



^ Since the surface brightness decreases rapidly with radius, the signal in the outer regions is dominated by low photon counts s and 
high background; consequently, the measurements are not as robust in the outskirts of clusters. 
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Figure 2. Creating a stable temperature profile. The black shows the spherically averaged temperature profile obtained for the gas in 
hydrostatic equilibrium. The red shows the temperature profile after an exponential drop-off has been added to the outer region. The 
dashed line shows the core radius, Tc = Rvir/3c for the gas and corresponds to the locatio n of the dip. Th e temperature drops from 
0.7T vir at 0.2Rvir to 0.4Tvir at 0.6Rvir and is consistent with observations of galaxy clusters (jVikhlinin et al.lboQS : Leccardi Molendi 
l2008h . 



Table 1: This table lists the parameters used to make the 9 isolated galaxy models. A 
constant physical softening of 1.14 kpc was used for both the gas and dark matter particles 
in all simulations. The number of particles used to model the two smaller galaxies was 
chosen such that the mass of each particle (dark matter and gas) was the same as the mass 
of the corresponding particle for the primary galaxy. 



Gas Content 


Galaxy Type 


Cone. 


Vcirc 


Ftvir 


Mvir 


Mdm 




Gas Mass Res.§ 








[km/s] 


[kpc] 


[10^° M0] 


[10^° M0] 


[10^] 


[10® M0] 




Primary 


10.0 


160.0 


228.6 


136.1 


150.0 


5.0 


1.5 


1% gas 


Onethird 


16.0 


110.0 


158.5 


45.3 


48.0 


1.6 


1.5 




Onetenth 


25.0 


74.0 


106.1 


13.6 


14.5 


0.48 


1.5 




Primary 


10.0 


160.0 


228.6 


136.1 


150.0 


5.0 


15.0 


10% gas 


Onethird 


16.0 


110.0 


158.5 


45.3 


48.0 


1.6 


15.0 




Onetenth 


25.0 


74.0 


106.1 


13.6 


14.5 


0.48 


15.0 




Primary 


10.0 


160.0 


228.6 


136.1 


150.0 


5.0 


27.0 


18% gas 


Onethird 


16.0 


110.0 


158.5 


45.3 


48.0 


1.6 


27.0 




Onetenth 


25.0 


74.0 


106.1 


13.6 


14.5 


0.48 


27.0 



^ The number of gas particles is the same as the number of halo particles in each galaxy. 
§ We used a fixed 50 SPH neighbours for all the simulations. 
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Figure 3. Evolution of gas density and temperature for a 10% gas fraction in our fiducial galaxy during initial condition construction. 
Left. The evolution for both the gas and dark matter is minimal; however, we see some outflowing features in the dark matter density. 
This results in a slightly larger galaxy at the end of 1 Gyr. The analytical input proflles shown in magenta (for dark matter) and cyan 
(for gas). The vertical dashed line is the softening used in the simulation. Right: The virial temperature is noted in the bottom-right. As 
the system evolves, the temperature proflle drops off exponentially at roughly Rvir with a 'knee' around ~ O.SRvir- 



Table 2: This table lists some of the physical quantities for the galaxies after they were 
evolved with zero artificial viscosity for 1 Gyr. The mergers are generated from two galaxies 
in this table. 



Gas Content 


Galaxy Type 


Unbound gas 


Lx 


Rvir 




Hot gas 


Virial of Clausius 






[%] 


[10^° erg/s] 


[kpc] 


[10^° M0] 


[%] 








Primary 


0.36 


0.10 


224.3 


128.9 


70.0 


1.11 


1% gas 


Onethird 


0.50 


0.03 


154.5 


42.1 


46.0 


1.08 




Onetenth 


0.58 


0.01 


103.9 


12.8 


21.0 


1.07 




Primary 


0.62 


8.87 


223.4 


126.5 


79.0 


1.12 


10% gas 


Onethird 


0.43 


3.17 


153.9 


41.6 


48.0 


1.08 




Onetenth 


0.56 


0.88 


103.7 


12.8 


24.0 


1.06 




Primary 


0.68 


26.5 


222.9 


126.3 


77.0 


1.12 


18% gas 


Onethird 


0.42 


10.9 


153.6 


41.3 


50.0 


1.08 




Onetenth 


0.50 


2.75 


103.3 


12.6 


25.0 


1.06 



2.3 Dependence on Artificial Viscosity 

Since the AVP determines the strength of the shocks produced in the simulation, we test the effect varying the AVP. We 
evolve the primary galaxy with a 10% gas fraction in isolation, with varying AVP from 10 ~^ — 1.0. Note that the recommended 
range of AVP is 0.1 — 1.0 for numerical simulations. We find that the difference in the final temperature profile is negligible 
between an AVP of 0.1 and 1.0; similarly for AVP's of 10~^ and 0. However, the average temperature is 20 — 30% larger in 
the higher AVP case. Fig. |4] shows the evolution of density, temperature and Lx from these set of simulations. In keeping 
with the recommended value for AVP, we chose an AVP= 0.5 to run all of our merger simulations. We also ran an equal-mass 
merger with 10% gas fraction and 0.1 Rvir impact parameter with three different AVP's of 0.1, 0.5 and 1.0. The results 
are indistinguishable between the three, showing that we are in a convergent regime for this physically reasonable value of 
AVP= 0.5. 



2.4 Merger Orbits 



Cosmological simulations of large-scale structure formation shows t hat ^ 40% of the mergers occur with an eccentricity 
e = 1.0 ± 0.1 and ^ 85% have an impact parameter b > 0.1 Rvir (i Khochfar &: Burkert 20061 : BensonI 2005). More than 
50% of the mergers occur with b > 0.5 Rvir- Keeping these larger impact parameters in mind, we create three orbits with 
b = 0.01, 0.1 and 0.5 Rvir; - two of our three orbits probe the larger impact parameters. We chose elliptical orbits with 
e ^ 0.05, 0.2 & 0.45. We also have a suite of equal-mass hyperbolic mergers. Our simulations contain a set of 27 elliptical 
and 9 hyperbolic orbits, for a total of 36 simulations. Table [3] lists the eccentricity, total energy and the initial velocities for 
each merger simulation. 
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Figure 4. Top Left: The gas (lower curve) and dark matter (upper curve) density profile with the analytical f3 and the NFW profile 
overlaid as solid lines, for the primary galaxy with gas fraction of 10%. The galaxy is represented by 5 x 10^ dark matter and gas particles 
each and is evolved for 1 Gyr with an AVP= 0.1. The fiattening of the dark matter profile occurs in the inner regions of the halo and 
corresponds to the softening - represented by the vertical dashed line. Top Right: The same plot showing the gas density profile at 1 Gyr 
for different AVP's. Bottom Left: The evolution of the temperature profile for different AVP's. The temperature increases by almost a 
factor of 2 for AVP > 0.01. Bottom Right: The evolution of the total X-ray emission via thermal bremsstrahlung from gas particles with 
T > lO^-^K and p < 0.01 Mq/pc^. 
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Table 3: The list of all the merger simulations performed. The first three 1:1 mergers for 
each gas fraction are the hyperbolic encounters, while the rest are all bound elliptical orbits 
with a fixed orbital energy for a given merger ratio. 



Gas 


IVIerger Type 




Eccentricity 


Orbital Energy 


Primary Vel. 


Secondary Vel. 


— 




[kpc] 


— 


[10^6 erg] 


[km/s] 


[km/s] 




1 : 1 


2.3 


2.99 


2864.6 


137.8 


137.8 




1:1 


22.8 


2.95 


2864.6 


137.8 


137.8 




1:1 


114.3 


2.78 


2864.6 


137.8 


137.8 


CO 


1:1 


2.3 


0.06 


-1432.3 


68.9 


68.9 


03 

bJO 


1: 1 


22.8 


0.18 


-1432.3 


68.9 


68.9 


1:1 


114.3 


0.41 


-1432.3 


68.9 


68.9 




3:1 


1.9 


0.06 


-542.0 


29.5 


92.2 


3:1 


19.4 


0.20 


-542.0 


29.5 


92.2 




3:1 


96.8 


0.44 


-542.0 


29.5 


92.2 




10:1 


1.7 


0.07 


-188.5 


10.5 


108.9 




10:1 


16.7 


0.21 


-188.5 


10.5 


108.9 




10:1 


83.7 


0.48 


-188.5 


10.5 


108.9 




1:1 


2.3 


2.99 


3397.9 


143.9 


143.9 




1:1 


22.8 


2.95 


3397.9 


143.9 


143.9 


oj 


1:1 


114.3 


2.78 


3397.9 


143.9 


143.9 


1:1 


2.3 


0.05 


-1698.9 


71.9 


71.9 


bJO 


1:1 


22.8 


0.18 


-1698.9 


71.9 


71.9 


1:1 


114.3 


0.41 


-1698.9 


71.9 


71.9 




3:1 


1.9 


0.06 


-642.9 


30.8 


96.3 


O 


3:1 


19.4 


0.20 


-642.9 


30.8 


96.3 




3:1 


96.8 


0.44 


-642.9 


30.8 


96.3 




10:1 


1.7 


0.07 


-223.6 


10.9 


113.6 




10:1 


16.7 


0.21 


-223.6 


10.9 


113.6 




10:1 


83.7 


0.48 


-223.6 


10.9 


113.6 



GO 



1:1 


2.3 


2.99 


3910.1 


149.0 


149.0 


1:1 


22.8 


2.95 


3910.1 


149.0 


149.0 


1:1 


114.3 


2.78 


3910.1 


149.0 


149.0 


1:1 


2.3 


0.05 


-1955.1 


74.5 


74.5 


1:1 


22.8 


0.18 


-1955.1 


74.5 


74.5 


1:1 


114.3 


0.41 


-1955.1 


74.5 


74.5 


3:1 


1.9 


0.06 


-739.8 


31.9 


99.7 


3:1 


19.4 


0.20 


-739.8 


31.9 


99.7 


3:1 


96.8 


0.44 


-739.8 


31.9 


99.7 


10:1 


1.7 


0.07 


-257.3 


11.3 


117.6 


10:1 


16.7 


0.21 


-257.3 


11.3 


117.6 


10:1 


83.7 


0.48 


-257.3 


11.3 


117.6 



3 ANALYSIS 

3.1 Calculating X-ray emission 

We assume that the primary source of X-rays in a galaxy merger is thermal bremsstrahlung from hot, diffuse gas of primordial 
composition. In real galaxies, X-rays are also produced by discrete sources like supernova remnants, compact binary objects 
as well as accretion events occurring around massive black holes. However, such emission is more likely to occur in the denser 
central regions of the galaxy and hence, also more likely to be obscured by the larger column density of the inter vening gas 



and d ust. The hot, diffuse gas should be relatively less affected by in situ absorption (jCox et al.ll2006 V Following ICox et al 



( 20061 ). we identify only those particles that have T > 10^'^ K and density ^ 0.01 M© pc ^ as X-ray emitting hot, diffuse gas. 



Throughout this paper, we will use the term hot gas to denote those gas particles that satisfy this criteria. Therefore, we 
estimate the X-ray luminosity as: 

Lx = 1.2 X 10-'^ (/imp)-' ^m,/>,7;'/'ergs-\ (4) 

i 

where rrip is the mass of a proton, mi,pi,Ti is the mass, density and temperature of a gas particle respectively and the 
summation is over only the hot gas particles 

The orbital energy in the gas is ^ 10^^ erg; much of this energy gets redistributed as the internal kinetic energy of the 
remnant. The average X-ray luminosity in all the simulations is ^ lO'^^ erg/s and it would take more than 10 Gyrs to radiate 
away only the orbital energy via bremsstrahlung. To ensure that we select only those gas particles that will remain hot, we 
compute the cooling times for the hot gas particles and reject the ones that have cooling time less than the dynamical time. 
As a further safeguard, we also reject the hot gas particles that have a cooling time less than the time between consecutive 
snapshots (~ 100 Myrs for all simulations). After these two step rejection process we recompute the X-ray emission from the 
hot gas and the shocked gas particles and find that the X-ray emission changes by less than 5% for the hot gas and less than 
1% for the shocked gas particles. Thus, our predictions regarding X-ray emission from the hot and the shocked gas are not 
susceptible to significant changes from radiative losses. 



^ Note that enriched gas at ~ 10^ K will cool primarily through metal line emission, implying that the X-ray emission could well be 
enhanced in actual colliding galaxies. 
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3.2 Shock Detection 



In our simulations, large-scale shocks develop in the gas. Since our simulations are adiabatic, shocks ar e the only mechanism 
that can increase the entropy of a gas particle. We take advantage of the conservative-entropy formulation ([Springel Hernauist 
of GADGET-2 ([Springell I2OO5I ) to identify the shocked particles based on their rate of entropy change. Each snapshot 



file in GADGET-2 gives the rate of change of entropy for every gas particle. We set a threshold for identifying shocked gas 
particles. We chose the threshold conservatively, i.e., we are likely to miss some shocked particles. The threshold, in code 
units, is is 10^°, 3 x 10^ and 2 x 10^ for 1%, 10% and 18% gas fractions respectively. We also locate every particle that was 
ever shocked by analysing all the snapshots and tracking the particle ids. This helps to determine the mechanism behind a 
particle's unbinding. 



3.3 Unbound particles 

At the beginning of the simulation, we identify the most bound 10% of the particles of each galaxy and track them to follow 
the centre of mass. The centre of mass velocity is subtracted to compute the total energy of individual particles. Particles 
that have a total energy greater than zero are considered unbound from the system. By defining the unbound material in 
this fashion, we avoid tagging particles that could be unbound from each galaxy but not unbound from the entire system. In 
the simulations where we have a merger remnant, we identify the unbound particles at the end of the simulation (when the 
system has settled down) and track those particles through the duration of the merger. 



4 RESULTS 

In this section we explore how different merger ratios, gas fractions and impact parameters affect the X-ray emission and the 
unbound gas. A priori, we expect to find an increasing trend for X-ray emission with increasing progenitor masses - this is due 
to the higher virial temperature and larger hot gas fraction in the more massive haloes (see Table |4|. In addition, a smaller 
impact parameter implies denser parts of the galaxy inter-penetrate; this should lead to stronger shocks and larger hot gas 
fractions. Therefore we expect the X-ray emission from both the regular hot gas component and the shocked component to 
increase with a decreasing impact parameter. A larger impulse is delivered during a close encounter with more massive galaxies; 
therefore, we expect an increasing trend in the unbound gas fractions with increasing galaxy mass and decreasing impact 
parameter. Thus, both Lx and the unbound gas fraction should behave in an analogous manner; we test these hypotheses in 
this section. 



4.1 X-ray emission 

Though the isolated galaxies have varying initial hot gas fractions, during the course of the merger more than 90% of the 
gas is in the hot phase. This is caused primarily by shock heating during the encounter. For example, the smallest secondary 
galaxies have ^ 40% of their gas in the hot phase in isolation (see Table |4|, but that fraction increases to more than 90% 
during the merger process. The total X-ray luminosity of the merger is determined by this hot gas fraction and can be seen 
from the correlated decrease in both the hot gas fraction and the X-ray luminosity (see Fig. [5] and Fig. [6}. 

For equal-mass mergers, the elliptical and hyperbolic orbits produce a peak Lx of ^ 10^° erg/s for the 1% gas fraction. 
This is from the first pericentre pass, which creates the strongest compressive forces. Since Lx is proportional to the square of 
the density, the shocks at this stage create the largest X-ray emission. However, note that even though we are talking about 
close passages, the average distance between the centres is greater than 20 kpc, greater than the typical disk size in spiral 
galaxies. Hence, the shocks will be observed at large distances from the centres. For equal-mass mergers with 10% gas, Lx 
peaks at ~ 8 X 10^^ erg/s during pericenter pass. The 18% gas fraction simulations produce Lx as high as - 2 X 10^^ erg/s. 

The simulations with the smallest impact parameters cause the denser parts to collide. As the impact parameter increases 
to 0.5 Rvir, the lower density regions of the galaxy inter-penetrate. Since our simulations are adiabatic, the maximum increase 
in the density can be a factor of 4; the strong dependence of Lx on the gas density results in a lower Lx • 

The equal-mass mergers create the strongest shocks in our entire suite of simulations. However, the surprising part 
is that the strongest temperature enhancements occur after the galaxy centres have passed by. At this point the shocked 
material encounters a lower density material upfront and accelerates forward. The material in the far side of the galaxy is 
still continuing in orbit with a velocity opposite to that of the accelerated shocked material. When these two fluids collide, 
the strongest shocks (measured by the rate of entropy change) are produced (see Fig. [7land[8|. This can also be seen from 
the temperature projections in Fig. [9] where the highest temperatures are created by the gas that is thrown forward after the 
pericenter passage, colliding with the gas in the far side of the galaxy. Consequently, the peak Lx from shocks occurs some 
time after the peak Lx • The equal-mass mergers with 1%, 10% and 18% gas have a peak Lx from shocks of ^ 10^^, ^ 10^^ 
and - 3 X 10^^ erg/s. 
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Figure 5. The evolution of the hot gas fraction for all the 36 simulations. The left column contains 1:1 mergers, the middle column 
contains 3:1 mergers and the right column contains the 10:1 mergers. The top row, middle and the bottom row show simulations with 
1%, 10% and 18% gas fractions respectively. The blue, red and the green shows the 0.01, 0.1 &; 0.5 Rvir impact parameters respectively. 
The dashed lines in the left column show the hyperbolic 1:1 mergers. We see that the peaks of the hot gas fractions corresponds to the 
pericenter passages. The most drastic evolution is seen in the equal-mass mergers where hot gas fraction evolves from ~ 80% to ~ 50% 
for 0.01 and 0.1 Rvir- For all the other simulations the final hot gas fraction reaches equilibrium at around 80%. 



The 3:1 mergers show a similar pattern; with Lx increasing with decreasing impact parameter and increasing gas fraction. 
The peak Lx for the 3:1 mergers with 1% gas is ^ 10^^ erg/s while the peak Lx from shocks is ~ 5 x 10^^ erg/s. The 10% 
gas fraction produces a peak Lx of ^ 10^^ erg/s and ~ 10^^ is from shocks. For the 18% gas fraction, the peak Lx and the 
peak Lx from shocks are ^ 10^^ and ^ 10^^ erg/s respectively. Though peak Lx from shocks is larger than 10"^^ erg/s it is 
still 3 orders of magnitude smaller than the total Lx signature from the merger event. 

The 10:1 mergers follow the trend and for the 1% gas fraction simulations the peak Lx and the peak Lx from shocks 
are ^ 2.0 x 10^^ and ^ 10^^ erg/s respectively. For the 10% simulations the peak Lx is 2 x 10^^ erg/s while the peak Lx 
from shocks is ^ 3 x 10^^ erg/s. For the 18% gas fractions, the peak Lx is ^ 5 x 10^^ erg/s while the peak Lx from shocks 
is - 2 X 10^^ erg/s. Clearly, the 10:1 mergers do not produce any appreciable X-ray signature from the shock-heating of the 
hot gas even with the largest gas fractions. See Table [6] for detailed data on X-ray production for all the simulations. 

From our simulations, we see that the expected trends do occur. The X-ray luminosity, both intrinsic and from shocks, 
increases with increasing halo mass and decreasing impact parameters. Larger initial velocities do generate larger Lx ; however, 
if the hot gas fraction drops too much, as it does in the case of some of the hyperbolic mergers, the X-ray emission can be 
reduced. In light of these trends, the best bet to detect the shock signature of hot halo gas during mergers would be to locate 
major mergers between Milky Way-type galaxies that are happening in the local universe. Such a system would produce 
significant X-rays from the hot halo and the shock-heated interface if the galaxies are to be observed in the initial stages 
of the merger. Spatial resolution will not be an issue since the signature is an extended source; the problem would be the 
identification of the underlying shock X-ray emission and separate it from the emission from the hot halo particles themselves. 
XMM-Newton, with its superior collecting area, would have been the instrument of choice; however, the inferior sensitivity of 
XMM-Newton in the soft X-ray band, where most of the shocked gas in our simulation radiates (the maximum temperature 
reached in our simulations is ^ 0.3 keV), Chandra is best suited to characterise the hot halo gas. Our threshold for the X-ray 
flux, 10""*^^ erg/s/cm^ is relevant to Chandra. 

So far we have discussed the production of X-rays by the merger, both from the hot gas in the halo and the shocks. We 
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Figure 6. The X-ray luminosity versus time for all simulations. The peak X-ray luminosity increases with increasing gas fraction oc p^^^ 
and occurs at the first pericentre pass. The peak X-ray luminosity also shows an increasing trend with decreasing impact parameter and 
decreasing secondary galaxy mass, correlated with the amount of shock heating that occurs. The line types and the colours are the same 
as in Fig [5] 



found that the X-ray emission depends on a variety of parameters like the merger ratio, the cumulative mass, the merger 
orbit and the gas fraction in the simulation. We attempt to combine the effects of all these factors to predict the peak X-ray 
luminosity from shocks as a function of the gravitational impulse from the merger. ICovington et al.l (|2008|) derive an empirical 
fitting formula for the impulse motivated by numerical simulations of galaxy mergers. The impulse, A^, is approximated by: 

^E(x ^ r ^ , T, (5) 

[^p +rp X Rvir,lJ 

where Mi is the primary galaxy mass, M2 is the mass of the secondary, and rp are the relative velocity of the galaxies and 
the pericenter distance as determined from the simulation. 

We also define a "merger strength", H, that is logically related to Lx • Lx depends on p^VT, and T oc ^^reiative (see 
Eqn. [1]). Therefore, Lx should depend on t'reiative; we use the velocity at perihelion, Vp determined from the simulation to 



capture this. Our final equation for the merger strength, AH, is then AE 



2 . 

gas ' 



and can be written as: 



AH oc 



[rl + rp X Rvir,i] 



(6) 



We then fit the peak Lx from shocks, Lpeak, as a linear function of AH. However, 
of 0.96, is obtained with a logarithmic fit. The fit for our set of 36 simulations is: 

-^^peak 



a better fit, with a correlation co-efficient 



logi 



0.88 X logio - L25. 



(7) 



]^Q38 erg/s 

This fit captures essential physics in a gas-rich galaxy merger and can be used to estimate the peak Lx from shocks in a merger 
(see Fig. [TO)) . Such a fit is quite useful in modelling the total Lx from shocks in semi-analytic studies of galaxy formation. 



4.2 Unbound Gas 



To constrain the amount of unbound material resulting from the mergers, we evolve the corresponding initial conditions in 
isolation for a similar time-scale. Every galaxy will lose some small amount of material from spurious two-body relaxation. 
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Figure 7. Temperature projection for the 18% gas fraction for a hyperbolic equal-mass merger with an impact parameter of 2.3 kpc. 
This run has the highest X-ray emission and largest unbound material. The X-ray contours are shown in purple, the unbound material 
with white contours and the portion of unbound material that is shocked with green contours. The shocked unbound only shows up in the 
4th image. The unbound material is identified at the end of the simulation and then traced through the entire run. Most of the unbound 
material comes from two distinct regions: the far lobes of the galaxy and the shocked front in between the galaxies. The projected X-ray 
emission also exceeds 10^^ erg/s in the shocked region with a total shocked X-ray emission exceeding 10^° erg/s. 

and by evolving the galaxies in isolation, we can estimate the magnitude of this effect. We deduct this amount of unbound 
material (column 3 in Tabled]) when we compute the final unbound fraction from the merger (see column 4 in Table [5]). 
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Table 4: Isolated galaxy evolution. We evolve isolated galaxies with AVP= 0.5 for ^ 10 
Gyr to ascertain the amount of unbound material from internal processes. This fraction 
is subtracted from the unbound material from the mergers, we can fix the actual amount 
freed during the merger itself. 



Gas Content 


Galaxy Type 


Unbound gas 


Lx 


Rvir 




Hot gas 






[%] 


[10^0 erg/s] 


[kpc] 


[10^° M0] 


[%] 




Primary ^ 


0.43 


0.09 


214.8 


113.2 


93.2 


1% gas 


Onethird 


0.69 


0.03 


149.1 


37.8 


73.4 




Onetenth 


0.84 


0.01 


100.2 


11.5 


41.3 




Primary 


0.60 


8.46 


212.9 


110.2 


89.8 


10% gas 


Onethird 


0.61 


3.38 


149.1 


37.7 


75.9 




Onetenth 


0.83 


0.94 


100.7 


11.7 


45.3 




Primary 


0.72 


27.5 


213.9 


112.0 


95.2 


18% gas 


Onethird 


0.60 


11.3 


148.6 


37.6 


79.0 




Onetenth 


0.76 


3.05 


100.7 


11.6 


48.4 



We find that most of the unbound material originates at or near the virial radius of the respective galaxies. As the merger 
progresses, this material escapes to beyond 3 Rvir even before the second passage. By the end of the simulation, this unbound 
gas is far 1 Mpc) from th e remnant (see Fig. [7|, thereby alleviating the problem of transporting gas to inter-galactic 
distances ( Ferrara et al. 200C)l l . We find that ~ 9% of the gas unbinds, along with a similar fraction in the dark matter. Per 
merger, this translates to a gas mass ^ 3 x 10^° Mq from equal-mass and ^ 10^° Mq from 3:1 mergers flowing into the IGM 
(see Fig. While at 18 Gyrs our simulation lasts much longer than a Hubble time, nearly all the unboun d material flows 



outside 3 i?vir within the first 5-6 Gyrs. Considering that every galaxy has had a major merger in the past (IConselice et al 



see), this may be an excellent mechanism to enrich the IGM. Thus galaxy mergers play an important role not only in the 
evolution of the galaxies and hierarchical structure formation, but they may also infiuence both the mass and the metallicity 
content of the IGM. 

We find that the unbound material originates from two distinct regions - from the near and far lobes of the galaxy along 
their line of motion. These constitute the regions that are shocked during the course of the merger. The near lobes are shocked 
before pericenter passage, which can be seen in Fig. [T] These shocked gas particles then leave in directions perpendicular and 
away from the direction of the motion of the galaxies as they expand into the lower density regions. Their escape is facilitated 
by the lower density and pressure of the ambient medium above the shocked region and can be seen from the velocity vectors 
of the unbound particles pointing outwards in Fig. [H The gas particles from this region escape from the galaxies the earliest. 

To understand why the gas escapes from the far regions of the galaxies, we must look into the shock process itself. As the 
galaxies pass through pericenter, the density in front of the shocked material suddenly drops rapidly. This causes the shocked 
material to accelerate, and it is fiung forward. However, the material at the far side of the galaxy is still moving with the 
bulk fiow. This naturally results in the largest relative velocity and hence, the strongest shocks in the simulation. The original 
shocked gas interacts with the gas from the far side, and the gas from the far side receives a large impulse, which causes the 
gas to fiow outward instead of continuing with the bulk motion of its host galaxy. This material can be seen as the brightest 
patches in the temperature images (see Fig. [8] and [9]). In unequal-mass mergers, the impulse affects the lower mass galaxies 
more, and an asymmetry is introduced whereby the strongly- shocked unbound material comes primarily from the smaller 
galaxy. This can be seen in the images of the temperature projections for the unequal-mass mergers where the brightest 
patches are towards the side of the smaller galaxies (see Fig. [9] for the 3:1 case). Although most of the unbound material 
originates from a radius greater than 0.5 Rvir in each galaxy and therefore should have a systematically lower metallicity, the 
gas can get enriched during the merger process itself from starburst- driven winds. 



t at 15 Gyrs. 
* at 10 Gyrs. 
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Figure 8. Temperature projections for the 18% gas fraction for a hyperbolic equal-mass merger with an impact parameter of 2.3 kpc. 
This run has the highest X-ray emission and largest unbound material. The velocities of a fraction of randomly chosen unbound particles 
is plotted; blue represents unshocked unbound material while white shows the shocked unbound material. We see that the unbound 
material escapes preferentially along the shocked front away from the plane of the merger since this provides the least resistance in terms 
of gas pressure. Strong shocks are generated once the centres have passed and material is ejected that interacts with the bound material 
at the other end of the galaxy that is yet to make a close pass (and hence has a velocity vector opposite to that of the unbound material). 
This is seen in the right image in the middle row. All of the unbound material leaves the vicinity of the galaxy much before they actually 
merge 6 Gyrs). 
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Figure 9. Projected temperature and Lx contours for a 3:1 merger with 18% gas fraction and an impact parameter of 1.9 kpc. A 
prominent bow shock develops as the secondary galaxy passes through the halo of the primary. Even at 2.2 Gyr, when the centres are 
separated by ~ 100 kpc, the strong shock persists and can be easily seen in the temperature map. Once the secondary material passes 
by the primary centre, a lot of material is ejected that gets shock-heated by the secondary galaxy as it moves inwards along the orbit. 
This can be seen in the patch of very bright material in the left-bottom image. The final image at 18.8 Gyr shows the merger remnant 
completely relaxed with spherical Lx contours. The minimum value of the Lx contour is lO'^^ '^ erg/s and increases in factors of 10°'^ 
erg/s. 
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100 10000 
Merger Strength 

Figure 10. The fit for the peak Lx from shocks as a function of the the merger strength, AH (see Eqn. [6]). Each simulation is plotted 
as a triangle while the fit (see Eqn. (T)) is plotted with a red line. The correlation for this fit is 0.96 and implies that Eqn. [71 is a good 
linear fit. 



18 Manodeep Sinha, Kelly Holley-Bockelmann 



Table 5: Unbound gas from mergers. The first three equal-mass mergers for each gas fraction 
are the hyperbolic encounters while the rest are all bound elliptical orbits with a fixed 
amount of energy in the orbit for a given merger ratio. The unbound fraction (Column 4) 
is obtained by deducting the unbound gas fractions when evolved in isolated (see Table |4t 
from the fraction determined at the end of the simulation. Column 5 shows the fraction of 
all the gas particles that are shocked during the simulation. The unbound gas particles are 
tracked through the snapshots to find where they become unbound or shocked. The fraction 
of unbound gas particles that are shocked during the simulation is shown in Column 6. 



Gas 


Merger Type 


b 


Unb. frac. 


Shocked (all) 


Shocked (unb) 


- 


- 


[kpc] 


[%] 


[%] 


[%] 




l:lh 


2.3 


19.2 


13.1 


49.6 




l:lh 


22.8 


18.5 


11.1 


46.7 




l:lh 


114.3 


11.2 


5.9 


44.1 




1:1 


2.3 


9.4 


14.1 


84.2 


03 
bJO 


1:1 


22.8 


8.4 


13.8 


85.8 


1:1 


114.3 


5.9 


11.9 


89.5 




3:1 


1.9 


6.7 


9.4 


69.6 


3:1 


19.4 


5.7 


9.2 


73.1 




3:1 


96.8 


4.1 


8.0 


80.0 




10:1 


1.7 


1.2 


2.6 


60.9 




10:1 


16.7 


1.1 


2.5 


62.6 




10:1 


83.7 


1.2 


4.0 


66.0 




l:lh 


2.3 


23.8 


22.5 


53.6 




l:lh 


22.8 


20.0 


13.4 


53.4 




l:lh 


114.3 


16.2 


8.4 


44.8 


1:1 


2.3 


9.7 


14.6 


86.9 


bJO 


1:1 


22.8 


8.9 


14.1 


86.8 


1:1 


114.3 


6.4 


12.1 
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36.4 
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5.3 
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It is clear that the unbound material is related to the extensive shocks in our simulations. Here, we attempt to disentangle 
the effects of the dynamics and the hydrodynamics of the gas as it unbinds. We backtrack the unbound material at the end 
of the simulation through the snapshots to identify when the gas particle first unbinds. At the same time, we note when 
and if the gas particle is shocked. Table [5] shows the results of this process. We find that for equal-mass mergers, galaxies 
on hyperbolic orbits systematically show that a smaller fraction of the unbound material originates from shocks, compared 
to galaxies on elliptical orbits. This may be because galaxies on hyperbolic orbits only undergo one pericenter pass, while 
galaxies on elliptical orbits experience multiple pericenter passes (and shocks). Elliptical orbits show a strong correlation 
between unbound and shocked particles. For instance, for equal- mass mergers, ^ 85% of all the unbound particles were 
shocked during the simulation. By tracking the particles back, we find that the overwhelming majority of gas was unbound 
right after being shocked. 

In a similar spirit as the previous estimati on of the peak Lx from shocks, we quantify the effect of the orbit on the 



unbound material. We follow the formulation of ICovington et al.l (|2008l l and use the impulse formulated in Eqn. [5] A linear 



regression analysis on /unb and and obtain a linear relationship of the form: 

/unb = 3.14 xlogioA^- 0.16. (8) 

The correlation coefficient for this fit is 0.95 (see Fig. I12[) . implying that the two quantities are well described by a linear 
function. Because we started out on a physically motivated basis by equating the impulse with the unbound fraction, this 
high correlation points to an underlying relation between the impulse and the unbound fraction. Since most of the hyperboUc 
mergers do not leave a remnant and the 10:1 mergers have only 1% unbound gas,, we do not use those simulations to fit this 
equation. Thus, out of our total 36 simulations, we have used only the 18 1:1 and 3:1 with elliptical orbits to obtain Eqn. [51 
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b= 0.01 0.1 0.5 




Time [Gyrs] 



Figure 11. Evolution of the unbound gas fraction for all simulations. Equal-mass mergers with the 0.01 &; 0.1 Rvir has ~ 20% unbound 
gas fraction from the final merger remnant. The 3:1 mergers release about 5-7% gas while the 10:1 mergers unbind the least amount of 
gas at about 0.5%. The irregular nature of the unbound gas fraction for the 3:1 and the 10:1 mergers is caused by the multiple passages 
taken by the secondary galaxy. 
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Impulse 

Figure 12. The linear fit (red) for the unbound fraction and the impulse, AE (see Eqn. [5]). We only use the 18 simulations with the 
1:1 and 3:1 bound orbits since they leave a merger remnant in all cases and have unbound fractions > 3%. We will explore this in an 
accompanying paper to find the unbound gas fraction throughout the history of structure formation. 



Table 6: Characteristics of the merger simulation and the final remnant for all the sim- 
ulations. Lack of an entry implies that the merger was not completed within a Hubble 
time. 



Gas Content 


Peak Lx^ 


Peak Shock Lx* 


Unb. gas mass 


Total unb. mass 


Hot Gas 




Gas within R-vir 


DM within Rvir 


Remnant Lx 




[10^*° ere /si 


r-|Q38 ere/sl 


[10-'-° Mol 


flO-'-^ Mol 

|_ X W -LVi J 


[XOJ 


[kpc] 


[XOJ 


[XOJ 


[10'*° ere /si 

[XU CI OJ 




0.82 


0.25 


0.57 


44.9 


51.0 








0.03 




0.54 


0.09 


0.56 


17.5 


55.2 


- 


- 


- 


0.04 




0.27 


0.005 


0.35 


9.3 


71.8 


- 


- 


- 


0.1 


m 


1.17 


0.19 


0.29 


2.03 


80.2 


248.1 


58.1 


58.4 


0.13 




0.81 


0.15 


0.27 


1.95 


82.1 


247.6 


59.1 


57.8 


0.11 


hn 


0.42 


0.03 


0.27 


2.6 


84.8 


246.3 


55.2 


56.7 


0.09 




0.50 


0.08 


0.14 


3.6 


84.2 


224.4 


62.9 


65.1 


0.07 


0.40 


0.04 


0.12 


2.5 


84.6 


224.4 


62.0 


65.2 


0.07 




0.20 


0.009 


0.09 


2.9 


83.6 


222.1 


56.8 


63.0 


0.07 




0.23 


0.004 


0.03 


2.0 


90.7 


218.9 


68.4 


72.8 


0.07 




0.19 


0.008 


0.02 


1.4 


90.3 


218.4 


67.4 
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0.07 




0.13 
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0.03 


1.6 
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39.1 


3.0 
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24.5 


52.9 








3.3 




23.1 


0.5 


5.1 


12.9 


59.3 


- 


- 


- 
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79.3 


18.7 


3.05 


4.36 


79.4 


251.3 


53.8 


60.4 


12.5 


bJO 


71.9 


10.2 


2.84 


4.16 
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248.6 


53.3 


58.3 


9.5 


34.5 


1.7 
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3.86 


83.3 
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57.1 


8.1 




44.0 


7.6 


1.4 


4.2 


82.2 
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56.0 


64.9 


6.2 


O 
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3.4 


1.2 


3.3 
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55.3 
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223.9 


52.2 


64.6 
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55.5 
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) gas 
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6.4 


14.7 
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27.6 
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30.7 


5.0 


6.1 


81.9 


250.4 


58.2 


59.7 


44.0 
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4.7 


5.8 
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57.8 


60.6 


35.0 
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7.1 


3.4 


4.9 


85.9 


248.6 


54.0 


58.3 


28.2 




193.5 


20.8 


2.3 


4.7 


85.6 


224.8 


59.5 


65.3 


22.2 


00 


108.4 


14.2 


2.1 


3.9 


86.2 


224.4 


58.5 


65.2 


23.2 




61.8 


1.8 


1.6 


3.9 


86.7 


223.9 


54.6 


64.6 


19.6 




84.5 


3.5 


0.39 


1.9 


91.6 


218.4 


62.6 


72.2 


18.4 




54.8 


1.7 


0.38 


1.6 


90.8 


216.6 


60.7 


70.4 


19.3 




35.8 


0.3 


0.45 


2.8 


89.1 
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21.4 



Lx is calculated by using only the particles that have Tgas > lO^ "^ K and p < 0.01 Mq/pc ^. 

Shock Lx is obtained by adding the Lx for the particles that exceed a threshold for The threshold, in code units, is 10^°, 3 x 10^ and 2 x 10^ for 1%, 10% and 18% gas 
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5 DISCUSSION 

In our simulations, we found that the hot halo gas shock-heats to temperatures ^ 10^"^ K in equal-mass mergers and ^ 10^ 
K in unequal-mass mergers. This is reflected in a strong temperature jump in the regions between the two colliding galaxies, 
well before pericenter passage, and persists until the first pass. The strongest shocks, and correspondingly the largest X-ray 
luminosities due to shocked gas, are created after the pericenter passage. For gas fractions greater than 10%, this Lx from 
shocks remains above observable thresholds of ^ 10^^ erg/s for at least a period of ~ 300 Myrs, and thus should be detectable 
in ongoing mergers with Chandra and XMM-Newton. Compared to the merger timescale, an observable X-ray shock is ^ 20 
times shorter. During the merger, the hot gas itself also radiates and can have an X-ray luminosity of ~ 10^^ erg/s for galaxies 
with a gas fraction as low as 10%. Herein lies the difficulty of detecting the shocked gas: the shocked gas is 1000 times less 
luminous than the non-shocked gas. To observe such X-ray shocks from mergers, one would have to look for a signal 3 orders 
of magnitude smaller than the global X-ray level and the effect only lasts for a timescale that is an order of magnitude 
shorter than the merger. Merger-induced star formation and AGN activity adds to the confusion by contributing X-rays from 
accretion disks. 

The redeeming feature of our predictions comes from the fact that the X-ray emission occurs even when the galaxy 
separation is large and occurs in a region where the X-ray background is low. Therefore, the ideal procedure to validate our 
predictions would be to observe close pairs of galaxies with an undisturbed morphology (indicating that the first pass has 
not occurred), with current space telescopes like Chandra and XMM-Newton. In such a case, there will be no contamination 
from merger-associated starbursts and the Lx signature, if observed, could be uniquely ascribed to shocked halo gas. Our 
simulations assume that the gas is always heated to the viria l temperature in g alaxies, irrespective of the galaxy mass; recent 



numerical simulations predict that that may not be the case (iKeres et al.ll2005[ ). Smaller mass galaxies, similar to the smallest 



logi 



0.88 X logio Aif- 1.25. (9) 



galaxy we simulated (with Vdrc = 74 km/s), are not expected to heat above 10^"^ K. In such a scenario, the predicted Lx 
from shocks will be much lower. We can use this prediction to constrain the theories of galaxy formation (both semi-analytic 
and from numerical simulations) against actual observations. 

A detailed analysis of the feasibility of the X-ray detection is beyond the scope of this paper; however, we will assume a 
conservative threshold of the X-ray flux as 10~^^ erg/s/cm^ and some representative distance to galaxy pairs of ^ 20 Mpc. 
This sets a detection threshold in the X-ray luminosity of 5 x 10^^ erg/s. We will find the hot halo itself exceeds this threshold 
for all the simulations but the shock X-ray luminosity does not always do so. We empirically model the peak Lx from shocks 
based on our set of 36 simulations as follows: 

IQ38 erg/s _ 

Note that we have modelled Lx completely by thermal bremsstrahlung, which is a fairly inefficient radiation mechanism at 
~ 10^"^ K. Using XSPEC, we find if that gas is at solar metallicity, the total emission will increase by more than an order 
of magnitude. Thus, our predictions about Lx should be considered a lower limit to the likely X-ray emission in actual 
observations of colliding galaxies for a given hot gas mass. 

Our simulations show that gas escapes from galaxies during the course of the merger, with most of the gas propelled 
outside 3 Rvir even before the final merger event. For equal-mass mergers, the unbound gas fraction ranges from 10-20% 
depending on the impact parameter and orbit type. The unequal-mass mergers result in 3-5% unbound gas, with similar 
fractions coming from each of the merging galaxies. All of this unbound gas, presu mably enriched with me tals, flows out 
into the IGM. Observations of the IGM at low redshift show the presence of metals ( Danforth Shiihl 20081 ): cosmological 



simulations can account for the metals since it only takes trace amounts of metals to leave a large-scale signature. However, 
about 50% of the baryons are actually pre sent in the IGM tFukugita Peebles 2004 ) in a temperature-density phase-space 



that is difficult to probe via observations ([Danforth Shul 



20081) . For about a decade, numerical simulations of structure 



formation have been struggling to obtain a prescription that allows such large amounts of gas to escape from galaxies without 
violating other observational constraints like the star formation rate. Our results show that mergers are an effective way 
to extract gravitational energy and cause significant mass-loss from galaxies. This process is similar, in principle, to the 
gravitational heating t hat has been incorporated for the first time in semi-analytic galaxy formation recipes in a recent paper 
( Khochfar Ostrikerl [2OO8 ) . Our simulations provide the empirical fits that can be incorporated into such semi- analytic 



simulations. Though the amount of gas has a strong dependence on the initial orbital energies, we see tha t each merger 



releases ^ 2 x 10 Mq of material into the IGM. Given that each galaxy suffers a major merger in its past (Cons elice et al 
the total a mount of unbound gas can becom e significant over the lifetime of the Universe. We calculate this in a 
companion paper ( Sinha Hollev-Bockelmann 20091 ) . 



We use those merger remnants with unbound fractions larger than 3% to obtain a prescription for the unbound gas 
fraction as a function of the merger impulse, AE, as follows: 

/unb = 3.14 X logio - 0.16. (10) 
Overall, both X-ray shocks and the unbound gas production are two facets of the same phenomena: the conversion of grav- 
itational potential energy into thermal and kinetic energy. Even though cooling, star formation and other baryonic physics 
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were not included in our simulations, the tremendous amount of available gravitational potential energy and the subsequent 
conversion of a fraction of it into thermal energy continuously heats the gas at a much faster rate than it can cool via any 
feas ible cooling mec hanism, at least for the temperatures and densities involved in our simulations. 

Cox et al showed that it is possible to repopulate the halo with gas (and metals) via mergers of disk galaxies; 

mass- loaded stellar winds from high mass stars, and AGN activity can also enrich the gas in the halo. We know that the 
metal content of the gas can strongly enhance the X-ray production, even overshadowing the contribution from thermal 
bremsstrahlung. For example, for solar metallicity gas at 10^ K, the emission for a MEKAL model is 2 orders of magnitude 
larger than that of thermal bremsstrahlung. As part of a future project, we will include metals in the merging galaxies to 
determine how this affects the X-ray properties. We expect that the boost from metal line emission will make the emission 
more detectable, and may allow X-ray shocks to be observed for galaxy pairs with larger separations. 

Our simulations were performed with vacuum boundary conditions, but in reality galaxies are continuously evolving 
objects within the background of an expanding Universe. Pristine gas is falling into the haloes, and the haloes themselves are 
growing more massive. Cosmological simulations are required to fully model the mergers within the context of an evolving 
background. As part of a future project, we will use cosmological simulations with dark matter and gas to explore the 
phenomena of the shocks in the mergers. 

An important implication of this work is that using our fitting formulae we can estimate the total emission in X-rays 
from shocks and the total amount of unbound g as as a function oi z over the history of the Univ erse. To do this, we will 
combine an extended Press- Schechter formalism ( Lacev Cole 19931 : Bower 1991 : Bond et al. 1991) with our empirical fits, 
in a semi-analytic approach. The semi-analytic estimate for the unbound gas fraction over the lifetime of the Universe can 
be compared against the current gas fractions of galaxi es as well as the baryonic content of the IGM. We will explore the 
evolution of the unbound gas fraction of the Universe in Sinha &: Hollev-BockelmannI ( 20091 ) . 
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